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Abstract 

Using a numerical simulation, we study the effects of ambipolar diffusion and 
ohmic diffusion on the magnetic field evolution in the interior of an isolated 
neutron star. We are interested in the behavior of the magnetic field on a 
long time scale, over which all Alfven and sound waves have been damped. 
We model the stellar interior as an electrically neutral plasma composed of 
neutrons, protons and electrons, which can interact with each other through 
collisions and electromagnetic forces. Weak interactions convert neutrons 
and charged particles into each other, erasing chemical imbalances. As a 
first step, we assume that the magnetic field points in one fixed Cartesian di- 
rection but can vary along an orthogonal direction. We start with a uniform- 
density background threaded by a homogeneous magnetic field and study 
the evolution of a magnetic perturbation as well as the density fluctuations 
it induces in the particles. We show that the system evolves through different 
quasi-equilibrium states and estimate the characteristic time scales on which 
these quasi-equilibria occur. 
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1 Introduction 



It is thought that the decay of the magnetic field in magnetars is the main 
source of its X-ray and 7-ray luminosity since these objects appear to radiate 
substantially more power than available from the energy rotational loss pp. 
The spontaneous decay of the magnetic field could occur through ambipolar 
diffusion, Hall drift and ohmic diffusion, which are non-ideal magnetohydro- 
dynamics processes that occur in thousands of years, compared to dynamical 
sound and Alfven time scales of milliseconds to seconds. Ambipolar diffusion 
promote a dissipative magnetic field advection, through the movement of a 
bulk of charged particles relative to the neutrons, the Hall drift is a non- 
dissipative advection of the magnetic field caused by the electrical current 
associated with it, and the magnetic field ohmic diffusion is a dissipative 
process caused by the electrical resistivity. The time scales of these process 
were estimated in Ref. 2. For classical pulsar magnetic field strengths, these 
were found to be longer than the lifetimes of these stars, and therefore un- 
likely to be observationally important. For the case of magnetars, it was 
found that magnetic field decay by these processes might be occurring [U [3] . 
However, a full understanding of these processes, their interactions, and their 
effectiveness in neutron stars is still lacking. The work of |2J was analytical, 
and therefore useful to identify general processes and identify the relevant 
time scales, but not to address the action of the identified processes in their 
full nonlinear development, and their interaction with each other. The full 
evolution of the magnetic field can only be addressed by numerical simu- 
lations. Recent ideal three-dimensional single-fluid magnetohydrodynamics 
simulations showed that, in a stably stratified star, a complicated, random, 
initial field generally evolves on a short, Alfven-like time scale to a relatively 
simple, roughly axisymmetric, large-scale configuration containing a toroidal 
and a poloidal component of comparable strength, both of which are required 
in order to stabilize each other pE]. It is interesting to study the effects of the 
non-ideal processes studied in Ref. 2 on the evolution of this configuration. 
As a first step in this direction, here we simulate the decay of a magnetic 
field, including the effects studied in Ref. 2, in a system where the magnetic 
field points in one Cartesian direction but varies only along an orthogonal 
direction. It is shown that the magnetic field evolves through different quasi- 
equilibrium states and we estimate the characteristic time scales where these 
quasi-equilibria occur. 
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2 General Physical Model 



We model the stellar interior as an electrically neutral and lightly ionized 
plasma composed of three moving particle species: neutrons (n), protons 
(p), and electrons (e). This is a slight extension of the model of [2], since 
these authors considered the neutrons to form a motionless background in 
diffusive equilibrium. We account for strong interactions between neutrons 
and protons by writing their chemical potentials as functions of both their 
number densities /ij(r, t) = ^({rijif, £)}), with i, j = n,p [5]. We consider the 
electrons as non-interacting particles, thus we write their chemical potential 
as a function of their own density, /i e (r, t) = /i e (n e (r, £)). The processes 
of [2] occur in thousands of years, compared to dynamical time scales of 
milliseconds to seconds, thus we use a slow-motion approximation in which 
we neglect the inertial terms in the equations of motion for the particles. 
The plasma species are described as three fluids coupled by collisions and 
electromagnetic forces, satisfying the equations of motion [21 E] : 

UiV/jii - rij— V-0 + riiqi(E H x B) 

& c 

Z)T«"i«j(^i-^)- C 1 ) 

In Eq. [U Vi is the mean velocity, fii/c 2 is the effective inertia of each 
species, if) is the gravitational potential, E and B are the electric and mag- 
netic fields, and the last term represents the drag forces due to elastic binary 
collisions, which damp the relative motions between the different species. 
The collisional coupling strengths are parametrized by the symmetric matrix 
jij, whose elements are inversely proportional to the characteristic collisional 
times between particles of species i against species j and generally depend 
on the local density and temperature. Neglecting the acceleration terms in 
Eq. [H leads to a set of differential equations that are first, rather than second 
order in time [21 16] . We introduce a slight modification in the equation of 
motion for the neutrons by introducing an artificial friction-like term of the 
form —n n av n . The basic role of this extra term is to allow a fluid element 
(containing all the species) to relax to the MHD equilibrium in a time scale 
that should be short compared with the time scales we are interested in, but 
long compared to the dynamical time-scales, 

= -n n V/i n - n n ^Vlf) - InjnnTLjiVn - Vj) - Tl n OiV n . (2) 
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We combine Eq. [T] for the charged particles, with the induction equation 
V x E = -(l/c)(dB/dt), to obtain 



OB 
~dt 



- = Vx 



[v n + v A + v H ) x B 



V x 



c 2 V x B 
Ana 



Ipnf^e 



l_ v / 7en/ip 

ec \ 7 en + 7 pn 



x Vi/j. 



(3) 



The first term on the right-hand side of Eqj3] shows that the magnetic 
field is transported by the sum of the neutron bulk velocity v n , the ambipo- 
lar diffusion velocity va and the Hall drift velocity vh- The second term 
represents the ohmic diffusion, where the electrical conductivity is 



a 



(4) 



Finally, the last two terms in right-hand side of Eqj3] represent battery 
terms. The ambipolar and Hall drift velocities are given by the equations: 



lpn{Vp ~ V n ) + 7 e „(V e - V n ) 

va = 

7en ~r 7pn 

Jjx B/c) - n c (V/x c + (/i c /c 2 )V^ 



(5) 



7en 7pw / -» -» \ 7en 7pn X _B 

= (f p — f e ) = — . (6) 

7en "I - 7pn 7en "I - 7pn 47ITi c 6 

In Eq. El we defined // c = /i e + whereas A/i = // c — //„ is the chemical 
imbalance. Weak interactions between the particles balance chemical poten- 
tial imbalances between the charged particles and neutrons. We define the 
difference between the rates, per unit volume, of the reactions p + e — > n + v e 
and n — > p + e + z/ e as Ar = T(p + e — > n + z/ e ) — Y{n — > p + e + u e ) = AA/i, 
where A generally depends on density and temperature [2]. The continuity 
equations for the particles are given by the equations: 

f/71 ■ — * 

^ + V.(rM*)=±A(A/i). (7) 
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In Eq. [3, the + sign corresponds to the neutrons and the — sign to the 
electrons and protons. The electrical neutrality condition is n e = n p = n c , 
and we define the baryon number density as Ub = n n + n c . If we add Eq. [7] 
for charged particles and neutrons, we obtain the conservation law for ub, 

+ V • [n n v n + n p v p ) = 0. (8) 

3 One-Dimensional Model 

In this section, we consider a one-dimensional model where the magnetic 
field in the neutron star points in one Cartesian direction but varies only 
along an orthogonal direction as B(f,t) = B z (x,t)z. We take this as a 
perturbation on a non-magnetized background in hydrostatic and chemi- 
cal equilibrium. The physical variables characterizing the background are 
labeled by the sub-index 0. We perturb this background by introducing 
a magnetic field, which causes small perturbations of the particle densi- 
ties. Thus, we write the particle densities as rii(x, t) = noi(x) + 8rii(x,t), 
where i = e,p,n, and 8rii(x, t) <C n 0i (x). The chemical potentials are 
fii{x,t) = Hoi{x) + 5fii(x,t), where 8fXi(x,t) = k iB 8n B {x,t) + k ic 8n c (x,t), 
kiB = (d[ii/dn B ) n0c , and k ic = (dfii/dn c ) noB , with k eB = 0. We use the 
Cowling approximation, thus, we neglect the perturbations of the gravita- 
tional potential respect to the background value, i.e., if) = if) . From Am- 
pere's law, j x = (c/47r)(V x B) x = = n c e(v px - v ex ), thus, v ex = v px = v cx . 
Using Eqs. EJ El EJ d and El we obtain the following dimensionless set of 
equations for the evolution of the magnetic field and the particles: 

dB z __ - d(v cx B z ) _d_ (_ dB z \ 
dt dx 10 dx \ dx J ' 

d8n B d 

— t^- = —qZ {n 0n Vnx + n 0c v cx ) , (10) 
d5n r d , 



'n Qc v cx ) - \n cs $ (k 0c 5n c + k 0B 8nB) , (11) 



dt dx 
where 



__ d ( k nB 8n B + k nc 5n c \ n 0c fi 0n d ( k cc 5n c + k pB 8n B ^ 
0n dx\ /Z 0n J n 0n dx\ JI 0n j 



i 9(b; 

n 0n dx 



(12) 
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dx 



n 0r , 



dx 



(13) 



In the above equations we have defined dimensionless variables as a = 
a/a s , with the sub-index s standing for a characteristic value of the corre- 
sponding dimensional variable. We set the characteristic values as: x s = L, 



B s 

Us 



B™ ax , t s = ax 2 /(n s k s ), v s = B 2 /(8nn s x s a), 5n s = B 2 /(8irn s k s 



n 



0/i 



r^max l, l,max 
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nB i 



'o '7s 



^max 
Jen 



X 
k, 



pn 



"+7en ax - On the other hand, we have defined r = c 2 /47rcr , k 0c = k cc +k nc , 
= k pc + k ec , ko B = k P B — k n B, n cs = n cs /n s . The dimensionless parame- 
ters T = B 2 /(87m cs n s k s ), I = a/(n s j s ), i? = a\ s x 2 Jn cs , and u = ar s /(n s k s ) 
define the importance of the different processes. 



4 Results 

4.1 Characteristic evolutionary time scales 

The Eqs. |9l [101 [U] describe the temporal evolution of the variables B z , 5n B , 
and Sn c respectively. The characteristic time scales on which this evolution 
occurs can be estimated from these equations. The two terms on the right 
hand side (RHS) of Eq. [9] give the characteristic time scale on which the 
magnetic field evolves due to the coupling with the particles and due to 
ohmic diffusion. On the other hand, the RHS of Eq. [10] gives the time scale 
on which a baryonic density perturbation evolves due to the motion of the 
particles. Finally, the two terms on the RHS of Eq. [TT]give the characteristic 
time scale on which a charged density perturbation evolves due to the motion 
of the particles and due to weak interactions between them. In the next 
paragraphs we give an order-of-magnitude estimative of each one of these 
characteristic time scales and describe the quasi-equilibria that are achieved 
along the evolution. 

4.1.1 Time scale to achieve magneto- hydrostatic quasi-equilibrium 

When a magnetic field is present in the system, the magnetic force pushes 
on the charged particles. At the first stages of the evolution, the collisional 
coupling between the charged particles and neutrons makes them move with 
about the same velocity, thus, v cx ~ v nx . Neglecting weak interactions, 
from Eqs. [TU1 HH we obtain 5tib / tiq B ~ Sn c /n 0c 1. The induced charged- 
particle and neutron pressure gradients tend to choke the magnetic force. 



6 



The magneto-hydrostatic quasi-equilibrium is reached when there is a close 
balance between these oposing forces (see RHS of Eq. [T2l . In the following, 
we estimate the time scale over which this close-balance occurs. Considering 
that n 0c <C n 0n , from Eq. [12] we obtain the associated velocity to a given 
perturbation 5ns as v nx ~ v cx ~ (1/aL) [k n B + (n 0c /noB)k n c] Sub- From the 
Eq. HOI we estimate the time-scale to produce this perturbation as 

L5n B x s Sn s ax 2 s 
tsn B ~ ~ = — r- = t s . (14) 

n BVnx n s v s n s k s 

We require this time scale to be the shortest, so we define it as our 
fundamental unit of time. Note that the non-physical parameter a in Eq 
controls how fast the magneto-hydrodystatic equilibrium is reached. 



4.1.2 Time scale for charged particles to reach diffusive quasi- 
equilibrium 

Now we suppose that the quasi-equilibrium discussed in the last section has 
been achieved (i.e., v cx = v nx + va x ~ va x )- Thus, only the charged-particles 
pressure gradients balances the magnetic force. The charged-particle diffusive 
quasi-equilibrium is reached when there is a close balance between these 
oposing forces (see RHS of Eq. [13]). From Eq. [T3j we obtain the associated 
velocity to a given perturbation 5n c as va x ~ ((k cc + k P B)8n c ) / {no n ^ cn L). 
Neglecting weak interactions, from Eqs. ITOl [TT1 and [141 we estimate the time 
scale to produce this perturbation, 

t 5nc ~ ^ = I » 1. (15) 
a e 

Note that the dimensionless parameter e is inversely proportional to the 
collisional strength parameter 7 cn between the charged and neutral particles, 
thus 7 S controls how fast the charged fluid reaches the diffusive equilibrium. 
On the other hand, we require e < 1, so ts nc ^>t s = 1. 



4.1.3 Time scale to achieve chemical equilibrium 

The characteristic time scale on which a chemical imbalance is erased by the 
weak interactions can be estimated from Eq. [TT] Neglecting the first term 
on the RHS of Eq. [TT], comparing the terms with 5n c , and using Eq. [T31 we 
obtain 

W ~ = i » i, (16) 

a\ s xj, § 
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We require d <C 1, so t^ ^> t s — 1. Note that, in Eq. [TH]the parameter 
A controls how fast the chemical equilibrium is achieved. 



4.1.4 Ohmic diffusion time scale 

The time scale on which the magnetic field evolves by ohmic diffusion can 
be estimated from Eq. [9j Neglecting the first term on the RHS of Eq. [9] and 
using Eq. [TU we get 

— Tigkg 1 . , 

tohmic ~ = = > 1. (17) 

ar s uj 

In Eq. [T71 the parameter uJ, which is directly proportional to the parame- 
ter r s , controls how fast the magnetic field decays by the electrical resistivity. 
We require uJ <C 1, so t ohmic ^>t s = 1. 



4.1.5 Time scale for the evolution of the magnetic field ignoring 
ohmic diffusion 

Here we neglect the ohmic diffusion term in Eq. [9] and assume the system has 
reached the magnetohydrostatic equilibrium (v nx ~ 0). In order to estimate 
the time scale for the evolution of the magnetic field, we have to distinguish 
between two cases. In the first case ts nc <C t^, and we want to describe the 
evolution of the magnetic field for times longer than t$ nc , thus, the charged 
particles will achieve diffusive equilibrium before chemical equilibrium is re- 
stored. The necessary ambipolar velocity to achieve a quasi-equilibrium state 
on which the charged density perturbations do not change can be estimated 
from Eq. [TT]as va x ~ (Sn c LX(k 0c + k 0B ))/n 0c . On the other hand, if we ne- 
glect the ohmic resistivity, from Eq. [9] we can estimate the time scale L/vax 
on which the magnetic field evolves as 



t B 



T T $ 

The opposite case corresponds to ts Uc ^> and we want to describe the 
evolution for times longer than t&^. In this case, the system achieves chemical 
equilibrium before the charged particles achieve diffusive equilibrium, thus 
during this time scale v a x 7^ 0. From Eq. [HI t B ~ L/va x and from Eq. [T31 we 
have, vax ~ B^/(n 0n n 0c 'y cn 87t~L). From Eq. [H]we can estimate the time scale 
for the evolution of the magnetic field as 

tB ~ -=r ~ (19) 
I I e 
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4.2 Evolution of a Linear Magnetic Field Perturbation 

In this section, we study the evolution of a linear magnetic field perturba- 
tion 8B Z . We write the magnetic field as: B z (x,t) = 1 + SB z (x,t), where 
5B z (x,t) <C 1. We start the simulation with an initial magnetic field per- 
turbation of the form SB z (x,0) = Acos(lnx), and a homogeneous parti- 
cle background: 5n c (x,0)/n 0c = 5n n (x,0)/n 0n = 0. For this perturbation, 
the characteristic length on which the magnetic force varies is L/(/7r), thus, 
t s = l/(/7r) 2 = 0.025. We define n = n 0c , k nc = k, g = nk cc , (3 = n 0c kpB- 
In the interior of a neutron star we have the typical values no„ = 10 38 cm -3 , 
n = 3.9 x 10" 2 , p = -0.0775, k = -2.962 and p = 1.508 [5J. In our analysis, 
we neglect the ohmic diffusion (ou = 0). We use the parameters: ZJ = 0, 
T = = 0.05, e = 0.005, A = -0.0001, and I = 2. For these parameters, 
tsn c = 5 and = 20. The characteristic time scale to reach the diffusive 
charged-fluid quasi-equilibrium can be estimated as the minimun between 
these two values. Since ts nc "C ^a^, the magnetic field significantly evolves 
on the time scale, t# ~ 1/(T t?) = 400. In Fig. 1, we show the system's 
evolution during the time scale to achieve magneto-hydrostatic equilibrium. 
We label the different instants of this evolution with progressive numbers, 
i.e., the initial condition is labeled with the number (1). We see from Fig. 
1 that at the end of this time scale, 5n c /no c ~ 5n n /no n , and the magnetic 
perturbation has not evolved significantly. In Fig. 2, we show the system's 
evolution during the time scale where the diffusive charged quasi-equilibrium 
is reached. The plots labeled with number 1 correspond to the final instant 
of the previous scale, i.e., the scale where the system achieved the magneto- 
hydrostatic equilibrium. We see in Fig. 2 that at the end of this time scale, 
5n n /n 0n ~ and 5n c /n 0c has significantly grown in order to balance the 
magnetic pressure. The magnetic perturbation has begun to evolve signif- 
icantly. It is also shown in Fig. 2, that at the end of this time scale, the 
separation between the lines is narrower, which shows that the charged dif- 
fusive quasi-equilibrium has been reached. In Fig. 3 we show the system's 
evolution during the time scale where the chemical equilibrium is restored 
and the magnetic perturbation decays. Again, the plots labeled with number 
1 correspond to the final instant of the previous scale, i.e., the scale where 
the system achieved the charged diffusive quasi-equilibrium. We see in Fig. 3 
that at the end of this time scale the density perturbations have been erased 
and the magnetic field perturbation has decayed. 
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5 Conclusions 



Using numerical simulations in one dimension, we studied the effects of some 
non-ideal MHD processes on the magnetic field evolution in neutron stars. 
We found that the system evolves through succesive quasi-equilibria, and we 
estimated the characteristic time scales on which these quasi-equilibria occur. 
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Figure 1: Magneto- Hydrostatic quasi-equilibrium. We have defined, b = 
SB Z , n b = 5n c /n Qc , n n = Sn n /n 0n . Dotted lines in Fig. lb. refer to the 
neutrons and full lines to charged particles. Label times are: (1): t — 0, (2): 
f = 0.013986, (3): f = 0.0279725. 
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Figure 2: Diffusive Charged Fluid Quasi-Equilibrium. We use the same 
conventions than Fig. 1. The label times are: (1,1*): * = 0.0279725, (2, 2*): 
t = 3.23304, (3, 3*): t = 8.57481, (4, 4*): t = 9.64316, (5, 5*): t = 10.7115. 
The upper asteristik in each number refers to neutrons. 
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Figure 3: Decay of the magnetic field and Chemical Equilibrium. We use 
the same conventions than Fig. 1 and Fig.2. The label times are: (1,1*): 
t = 10.7115, (2, 2*): * = 116.83, (3, 3*): * = 329.067. 
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